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Abstract. A problem of finding stationary states of open quantum systems is 
addressed. We focus our attention on a generic type of open system: a qubit coupled 
to its environment. We apply the theory of block operator matrices and find stationary 
states of two-level open quantum systems under certain conditions applied both on 
the qubit and the surrounding. 
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1. Introduction 

In quantum mechanics, tlie density operator p of a quantum system is called a stationary 
state if [H, p] = 0, where if is a given time-independent Hamiltonian of the system. 
Since p satisfies the Liouville-von Neumann equation idtp = [H,p], it is clear that 
stationary states are invariant with respect to the transformation p UtpU} , where 
Ut = exp{—iHt) is the time evolution operator. In other words, stationary states do not 
change during the time evolution. 

For low- dimensional closed systems, the stationary states can be obtained relatively 
easily [1,2]. It is a common situation, that a small quantum system is immersed in other, 
mostly large, system called the environment [6]. Such an open system does not evolve 
unitarily in time. An analysis of open quantum systems [3~5] is much more complicated 
as they are a stage of a variety of physical phenomena [7-9]. The famous decoherence 
process [10] may serve as an example. In open quantum systems character of potentially 
existing stationary states is not obvious. 

There are various physical problems related to the properties of open quantum 
systems, which has already been addressed and intensively discussed (see e.g., [11,12]). 
Nevertheless, the procedure of deriving stationary states is not one of them. The 
existence and properties of stationary states have significant importance in quantum 
information processing and quantum theory itself. One can pose natural questions: 

(i) do the stationary states exist for a given open system? 

(ii) what features of a given model are responsible for existence of such states? 

(iii) how such states can be constructed? 

The answers to the above questions are still incomplete. For example, it is known 
that the stationary states exist for completely positive (CP) [13] evolution of the open 
system, this fact follows directly from Schauder fixed point theorem [14]. However, this is 
only an existential result and so far there are no available methods to determine explicit 
form of stationary states. Furthermore, the very existence of the stationary states in 
general case is an open problem e.g., in the presence of initial system-environment 
correlations [15]. 

The purpose of this paper is to propose the method of calculating the stationary 
states in the case of two-dimensional open quantum system. The theory of block 
operator matrices [16-18] is adapted to achieve this goal. In particular, we use the 
Riccati operator equation [19] to solve the eigenproblem for the total Hamiltonian. It is 
shown how to derive the stationary states by using the solution of the equation. 

2. Block operator matrix approach 

We begin with a brief review of the block operator matrices approach to the problem 
of decoherence in the case of a single qubit [21-23]. Let H be the Hamiltonian of the 
total system. We will assume that it has the following form 
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H = iJQ ® Ie + Iq ® i^E + Hint, (1) 

where Hq and He represent the Hamihonian of the qubit and the environment, 
respectively, while Hint specifies the interaction between the systems. The Hamiltonian 
H acts on the Hilbert space "Htot = ® "He, where "He is the Hilbert space (possibly 
infinite-dimensional) related to the environment. Iq and Ie are the identity operators 
on and "He, respectively. 

Since the isomorphism C^^He — "He ©"He holds true, the Hamiltonian (1) admits 
the block operator matrix representation [18]: 



H 



V 



on Via) = {v{H+) n v{v^)) © {v{v) n v{h_)) 



(2) 



All the entries of (2) are operators acting on "He- Moreover, the diagonal entries, i.e., 
H± are self-adjoint. In this paper, we will focus on the case in which V is bounded, 
thus V'^ is bounded as well; however, no assumption on boundedness of is made. 
Under these circumstances we have I^(H) = V{H^) © V{H_), where domains V{H^) 
are assumed to be dense in "He- 

The generally accepted procedure to obtain the reduced time evolution of the open 
system, the so-called reduced dynamics, reads 



Pt = TTE[VMpo)Vl]=Ttipo). 



(3) 



Above, po specifies the state of the open system at t = 0. The map $ assigns to each 
initial state po a single state $(po) of the total system. The assignment map must be 
chosen properly so that Tt can be well-defined [24-26]. For instance, if no correlations 
between the systems are initially present, then $(po) = Po ® i^, for some initial state of 
the environment u. It is worth mentioning that, if the initial state cannot be factorized, 
the definition of $ is not accessible [27,28]. The unitary operator = exp(— iHt) 
describes the time evolution of the total system. 

The map TrE(-) : T('He © "He) — ^ M2(C) denotes the so-called partial trace: 



Trp 



Mn 


M12 




M21 


M22 





TrMn 
TrM2i 



TrMi2 
TrM22 



(4) 



Tr(-) refers to the usual trace operation on "He, T{1-Le ©"He) denotes the Banach space 
of trace classes operator with the trace norm: = Tr(Vv4At), whereas M2(C) is the 

Banach space of 2 x 2 complex matrices. Note, the partial trace is a linear operation 
transforming the block operator matrices (square brackets) to the ordinary matrices 
(round brackets). 
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3. Main results 

Fixed point theorems like Banach or Schauder indicate the existence of stationary states 
for a given evolution Tj. However, there is no general analytical procedure to obtain 
explicit form of such states. In this section we propose a method of deriving stationary 
states for two-level open quantum systems. The generalization to the higher-dimensions 
seems to be possible. However, we will not deal with this issue in this paper. We begin 
with some definitions. 

Definition 1. The density matrix p is said to be a stationary state if it is invariant 
with respect to reduced evolution, Tt{p) = p. 

Definition 2. Let X be an operator acting on the Hilbert space T-Le- The subset Tx of 
"He © "He defined as 

: e V{X) C -He^ (5) 

is said to be the graph of X. 

The graph of a linear and closed operator is a subset of the Hilbert space, which is 
a Hilbert space itself equipped with the inner product 

(^l|^2) = (^l|^2) + (0l|02), = 

{ip\(t)) is an inner product on "He- It is a known fact (see Lemma 5.3 in [20]) that the 
graph Tx is H— invariant, that is H (Fx H X'(H)) C Fx if and only if X is a bounded 
solution (with Ran(X|x){_f/+)) C V{H_)) of the Riccati equation: 

XVX + XH+- H_X -V^ = on V{H+). (7) 
Along with the equation above we introduce the dual Riccati equation, namely 

YV^Y + YH_- H+Y -V = on V{H_). (8) 

It is proved in [20] that Y = -X^ is a solution (with Ran(Xl"| C V{H+)) of (8) if 
and only if the orthogonal complement of Fx, i-c, the subspace 

: 1^) Gl)(Xt) cHe| (9) 

is H-invariant. It is straightforward to see that a bounded operator X solves (7) if and 
only if —X'^ is a solution of (8). Therefore, Fx as well H— invariant if and 

only if X is a bounded solution of (7). In other words. Fx is reducing subspace of H if 
and only if X is a bounded solution of (7). From considerations above follow also that 
Fx and F^ are U^— invariant. 




e F 



X 



1,2). 



(6) 




Stationary states 



5 



Definition 3. Elements from the graph and its orthogonal complement are denoted 
by \X^) and {X"^) , respectively. The Riccati states are defined as = Ti-Eig^p) CLnd 
= TrE(^'^), where := \X^){X^\ and := \X^){X^\. 

The vectors \X^) and are not normalized with respect to the norm induced 
by the inner product (6). Moreover, the states g^ and g^ are not factorisable {i.e., 
correlations occur), unless X\il)) ~ lip) and X^tp) ~ respectively. However, 
they are Uj(-)U|— invariant, which is obvious because the vectors \X^) and \X^) are 
Uf— invariant. As a consequence, the Riccati states and p^ are Tj— invariant, where 
the map Tf has been defined in (3). Therefore, the set of all the Riccati states is invariant 
under the time evolution. Nevertheless, the Riccati states are not the stationary states, 
in general. However, we show that the latter can be found among the Riccati states. 
To be specific, we will prove the following 

Theorem 1. Let X be a bounded solution of the Riccati equation (7). Then 

i) the Riccati state p^ is a stationary state if the vector is a eigenvector of the 

operator Z+ = H+ + VX : V{H+) He, 
a) the Riccati state p'^ is a stationary state if the vector 

operator Z_ = H_- V^X^ : V{H_) He- 



is a eigenvector of the 



Proof. Let Z^ 



Al 



for A G C and G V{H+). From (7) we obtain that 



+ H_X = XZ^, hence in view of (2) the last equality leads to H|X^) = A|X^). Thus, 
the vector state \X^) is the eigenvector of the total Hamiltonian with the corresponding 
eigenvalue A. Since H is self-adjoint we have A G M and in consequence Ut^^U 
which ultimately leads to Tt{p^) = p^. 

In a comparable manner we have H|X''^) = C,\X'^) for ,^ G M and 
that Z_\(f)) = ^10). Just as before Tt(^'^) = g'^, therefore Tt{p'^) = p^. 



G V{H_) so 
□ 



At this point, some remarks, regarding theorem given above, should be made. 

Remark 1. The question whether all stationary states are Riccati states or if it is 
possible that stationary states exist that are not Riccati states is still open. 

Remark 2. Since the space Fx is closed, we have the following decomposition T-itot = 
Tx®Fx- Thus the total Hamiltonian is similar to certain block diagonal operator matrix, 
S~"^HS = H(^, where 









z_ 



with V{Z±) = V{H^) and S 



Ie -Xt 
X Ie 



(10) 



This, implies that cr(H) = (j{Z+) U a{Z_ 
the eigenvalues of the Hamiltonian H. 



Therefore, the eigenvalues of Z± are exactly 



Proof. Let X be a bounded solution of (7). V is assumed to be bounded as well. From 
the definition of Z± we have T>{Z^) = V{H^), and thus 'C(H) = P(Hd). Since X solves 
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the Riccati equation (7), it is clear that HS = SHj. To prove H ~ we will show 
that S is invertible and is bounded. Indeed, S = I + X, where 



X 




X 







Since X^ = —X, the spectrum of X is a subset of the imaginary axis. In particular, 
— 1 ^ o"(X) thus ^ c"(S) and, hence, S has a bounded inverse. □ 



Remark 3. The stationary states p^, p"^ indicated in theorem 1 are given by 



A\ ^ ^^^^ 



and p^ = B 



(12) 



1 

where G V{H_^_) and |0) G V{H_) are the eigenvectors of Z_ and respectively. 
A = Tr(p^), B = Tr(p'^) are the normalization constants and (X)^ = {(p\X\(f). 

Proof. Since Tr|'(/')(0| = {(f)\ip), the equations (12) can be obtained directly from the 
definition (3) and the formula (4). □ 



4. Examples 

4.1. Spin-boson model 

In this subsection we will demonstrate an application of the presented method to a non- 
trivial example, namely, the paradigmatic spin-boson model [29,30]. Assume that the 
Hamiltonian of the qubit (spin-half) and its environment (boson) are in the following 
forms 



Hq = (3az + aax and He = uja'a, (13) 

respectively. For the sake of simplicity, we consider only the case of a single boson. The 
interaction between the systems reads 

Hint = ® {9*a + ga^) =a,^V. (14) 

In the above description, a^, cr^ are the standard Pauli matrices and a, /3 G M. The 
creation and annihilation a operators obey the canonical commutation relation (CCR) 
[a, a^] = Ie [1]. Parameters u > and (7 G C represent the energy of the boson and the 
coupling constant between the qubit and the environment, respectively. 

If a = (no energy exchange between the systems), the model can be solved, 
i.e., the reduced dynamics can be obtained exactly [31,32]. The solution describes the 
physical phenomena known as the pure decoherence or dephasing [33]. On the other 
hand, when a ^ the exact solution in not known. The objective is to estimate the 
stationary states for the latter case. 
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To proceed, we must clarify some technical aspects {e.g., domains of H^). Clearly, 
the operators a and a''' cannot both be bounded since the trace of their commutator 
does not vanish [31]. Therefore, the CCR holds only on some dense subspace V2 of He- 
Let Vi be the dense domain of both a and a\ on which they are mutually adjoint, that 
is, (a^)* = a, a* = a'^. At this point, it is not obvious that the sets Vi, V2, having desire 
properties, exist. The detailed construction can be found in [35] and the right choice is 
given by 

= I I^A) G He : X^n'=|(^|0n)|' < 00 I , k = l,2. (15) 

I 71=0 J 

On Vi the creation and annihilation operators can be defined explicitly as (see 
also [36-38]) 



= 5^V^(0n|0)|0n-l), at|0) = J2V^{(pn\mn+l), \4>) G V,. (16) 
n=l n=0 

{\(f>n)}'^=o is an orthonormal basis in T-Le- From (16) follows that a^cpn) = + l|0n+i) 
and a|0n) = \/n\(f)n-i), which in most books on quantum mechanics is a definition of 
the creation and annihilation operators. However, the operators defined in such a way 
are not closed; nevertheless, they are closable and their closures are given by (16). 

Since a is closed. He is a positive self-adjoint operator and V2 C Pi is a core of a 
(see e.g., Theorem 4.2.1 in [39]). Henceforward, we assume that the basis {|0n)}5^Lo is 
composed with the eigenvectors of He- In this case we have 

00 

HE\<P) = Y,^n{<P^\<P)\<Pn), |0)GI?2. (17) 

n=0 

By choosing a suitable coupling constant g, it is possible to make = He ± V self- 
adjoint on := VinV2 = T>2. To see this, let us first note that {g*a+ga'^)* D g*a+ga\ 
thus V is Hermitian (symmetric). Since He is self-adjoint and V Hermitian, it is 
sufficient to show that V is relatively bounded with respect to He {He bounded) and 
with HE-houYid less than one. The fundamental result of perturbation theory, known 
as the Kato-Rellich theorem [3 1] assure Hermiticity of H± in this case. 

Recall that S is A bounded if i) I? (A) C V{B) and ii) < a||A|0)||2+6|||0)f, 

for all 10) G 'D{A) and some nonnegative constants a, b. The infimum of all a for which 
a corresponding b exists such that the last inequality holds is called A— bound of B. 
Note, sometimes it is convenient to replace the condition ii) by the equivalent one: 

\\B\<j>)\\<a\\A\<l>)\\ + b\\m- 

It is not difficult to see that if two operators Bi, i = 1,2 are bounded with respect 

to the same operator A and their relative bound are less than bi and 62, respectively; 

then aiBi + 02-82 is also A bounded and its relative bounded is less than |ai|6i + \a2\b2 

(see Lemma 6.1 in [10]). In other words, the set of all A bounded operator form a linear 
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space. Therefore, to see that V is He bounded it is sufficient to prove that both a and 
are He bounded. To finish this note 



oo oo 

2 



5^n|(0„_i|0)p = ^(n + l)|(0„|0)r (18) 



n=l n=0 

oo oo 

n=0 n=0 



In comparable manner one can also verify that < Since 

if bound of both a and a'^ is less than one, the ifg— bound of V is also less than 
one for \z\ < 1/2. 

The block operator matrix representation of the spin-boson Hamiltonian reads 



H 



H+ a 
a H 



where H^ = He±V and ^'(H) =V2® V2] (19) 



the quantity a is understood as al-^- For the sake of simplicity we have set /3 = 0, the 
example remains non-trivial because [Hq ^ Ie, Hint] 7^ 0. The corresponding Riccati 
equation takes the form 

aX^ + XH+ - H_X - a = on V2. (20) 
In order to solve this equation we define an operator P as 

00 

P|^) = ^e^-(0„|v^)|0„), 1^) e He. (21) 

n=0 

Directly from (21) we have P = and P^ = Ie, thus P is both self-adjoint and unitary. 
Formally, P can be written as P = exp{i7rHE), however, unlike He, P is everywhere 
defined. The Hellinger-Toeplitz theorem guaranties that P is bounded, which can also 
be seen directly. Indeed, form unitarity we obtain \\Pip\\ = W^PW, for G He, hence 
||P|| = 1. P is, in fact, the bosonic parity operator [11]. We will show that X = P 
solves (20). Since P^ = Ie it is sufficient to show that 

PH - H+P = or equivalently PH^P = H^ on V2. (22) 

In order to prove (22) let us first note that Plip) G V2, for {tp) G V2, which means 
Ran(P|x)2) C V2. This follows from |e*™| = 1. Furthermore, PHeP = He and 
PVP = —V . The first equality is obvious, while the second one follows from PaP = —a 
and Pa'^P = —a^. And as a result we obtain (22). 

The stationary states in this example can be written as 



where r± = {ip^\P\tp^) (23) 
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and are the eigenvectors of = ± aP. Unfortunately, the solution of this 
eigenproblem is not known for a ^ 0. However, one can determine certain bounds on 
r± using properties of P and p±. First, r± are real numbers because P is self-adjoint. 
From the non-negativity of p± we obtain that r± G [—1, 1]. 

Equation (23) provides estimation of the stationary states of the spin-half immersed 
within the bosonic bath. This result has been obtained without any approximations. 
Of course, r± can be computed approximately with the use of known methods. It is 
important to stress that to obtain the exact reduced dynamics for the model in question 
one needs to resolve an eigenvalue problem for Z^. 

4-2. Commuting environment 

In the second example we consider the Hamiltonian in the following form 

li = aa^0lE + lQ^ Ho + a,^ Hi, a ^ 0. (24) 

Here we assume that linear, bounded operators Hq, Hi commute. In this example we 
impose restriction to the spectra of Hq and Hi, i.e. a{HQ), (y{Hi) are discrete and non- 
degenerated. The Hamiltonian (24) describes a qubit in contact with an environment 
and in the presence of the magnetic field B = Bcx, where B a. Examples of such 
systems occur in the literature, e.g. [42-44]. The block operator matrix representation 
of (24) is given by (19) with H^ = Hq±Hi and 'P(H) = "He- The corresponding Riccati 
equation reads (20). 

Using the fact that Hq and Hi commute, so they have a common set of eigenvectors, 
we write 

Ho\<i)n) = \n\(t>n) and Hi\(t)n) = ^n\(t>n) , (25) 

where A„ G a^Ho), ^„ G cr{Hi) and (0n|0m) = ^nm, for n, m G N. The Riccati equation 
has a positive and self-adjoint solution X = f{Hi), where the function / is given by 



/(a;) = ^/?±^L^, for xeaiHi). (26) 

a 

Unlike to the spin-boson model, in this case the eigenproblem for Z± can be readily 
solved. Indeed, we have 

Z±|0„) = [A„±e„/(en)]|0n). (27) 

According to the Remark 3 we obtain 

where C„ = (1 -|- |/(^n)P)~^- In this case there are no initial correlations between the 
systems since ~ |0n)- We wish to emphasize that there may exist other solutions 

of the Riccati equation. For an explicit example see Ref. [23]. 
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4-3. Sylvester equation 

In the case = 0, the Hamiltonian (2) is already in a block diagonal form. One can 
notice that the Riccati equation simplifies to the Sylvester equation: 



XH+ - H_X = 

There exist at least one solution, namely X = 
are given by the projections Pq = diag(0, 1) = 

4.4- All unbounded entries 



on V{H+). (29) 

: 0. The corresponding stationary states 
|0)(0| and Pi = diag(l,0) = 



In the last example we consider an interesting example in which all the entries of H are 
unbounded, but still the solution of the Riccati equation exists as a bounded operator. 
To see this, let us choose = Hq with domain V{Hq) and let us assume that V is 
self-adjoint with domain T>{V). Then, the Riccati equation 

XVX + XHo- HoX -V = on V{Ho)r\V{V) (30) 
has at least two bounded solutions, that is X± = ±Ie. The stationary states read 

P± = \ (^^ ' wJ^ere e± = ±1. (31) 

An example in which all the entries of a block operator matrix are bounded and the 
solution of the Riccati equation is unbounded has been provided in [16]. 



5. Summary. 

In this paper we have proposed the method of calculating the stationary states for two- 
level open quantum system. The theory of block operator matrices has been adapted 
to achieve this purpose. In particular, we have related the solution of the algebraic 
Riccati equation to stationary states. In the presented method, the stationary states are 
generated from the stationary states of the total system by tracing out the environment. 
Our investigation includes the case when the initial system-environment correlations 
occur. In fact, this case is embedded in the method since the eigenstates of the total 
Hamiltonian are entanglement. 

At the end, we want to stress that the method cannot be used when the total 
Hamiltonian is not known. Such a situation arises e.g., when the details of the interaction 
between the systems are not accessible. We hope that despite aforementioned weaknesses 
the results of the paper may serve as a starting point for further investigations. 
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